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摘 E: 基于 内 蒙古 45 个 气象 站 点 1980 一 2019 年 日 均 地 表 温 度数 据 .结合 中 国 第 一 代 全 球 陆 面 再 分 析 产 品 
(CRA) 数 据 以 及 NDVI 数 据 , 利 用 趋势 分 析 法 .相关 性 分 析 法 和 灰色 关联 度 , 对 内 蒙古 近 40 a 地 表 冻 融 指数 时 空 变 
化 特征 及 驱动 因素 进行 分 析 。 研 究 表明 :(D) SFI( 地 表 冻 结 指数 ) 年 均值 的 空间 分 布 特征 表现 出 自 西 南 向 东北 递增 
的 规律 ,STI( 地 表 融 化 指数 ) 则 反之 ,纬度 是 影响 地 表 冻 融 指数 空间 分 布 的 关键 因子 。 研 究 期 间 SFTI 和 STI 分 别 呈 现 


出 


(10a) FI 152.4 C- d° (10a); 与 季节 冻 土 区 相 比 ,多 年 冻 了 
近 40 a 土 壤 表层 含水 量 .降水 量  NDVI 呈 增加 趋势 , 雪 深 呈 减 少 趋势 ,但 年 际 变化 表现 出 不 同 的 空间 差异 性 ,多 年 


显著 下 降 和 上 升 趋势 ,多 年 变化 范围 分 别 为 956.1~1848.3 %.d 和 3717.6~4442.3 对 .d, 变 化 率 分 别 为 -156.4%.:d， 
上 区 的 冻 融 指数 对 气候 变 暖 的 响应 更 加 敏感 。(2) 研究 区 


冻 土 区 呈 暖 干 化 发 展 趋势 ,季节 冻 土 区 呈 暖 湿 化 发 展 趋势 。(3) 地 表 冻 融 指数 与 影响 因素 以 负 相关 关系 为 主 ,SFT 与 
影响 因素 在 多 年 冻 土 区 大 部 呈正 相关 关系 ,在 季节 冻 土 区 大 部 呈 负 相关 关系 ,STI 则 反之 。 内 蒙古 地 表 冻 融 指数 变 
化 受 影响 因素 共同 驱动 ,0.4m 土 壤 仿 水量 是 影响 SET 变化 的 主导 因素 ,NDVI 是 影响 SITI 变 化 的 主导 因素 。 研 究 结 


果 可 为 内 蒙古 冻 土 退化 . 农 牧 业 生产 等 提供 科学 的 参考 。 
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冻 土 一 般 是 指 温 度 在 0% 或 0% 以 下 ,并 含有 
冰 的 各 种 岩 土 和 土壤 ,包括 永久 冻 土 (PF) 和 季节 性 
冻 土 (SFG)。 我 国 是 世界 冻 土 大 区 ,多 年 冻 土 和 季节 
冻 土 分 别 占 国土 面积 的 22% 和 53.5%" IPCC 5 
次 评估 显示 ,全 球 地 表 温 度 在 1880—2012 年 期 间 升 
高 了 0.85 "2 。 全 球 变 暖 在 高 纬度 和 高 海拔 区 域 加 
剧 ,对 冻 土 的 形成 有 较 大 的 影响 ,尤其 是 在 我 国 北方 
高 纬度 ,温度 升 高 幅度 大 于 全 国平 均 升 温水 平 ? 3 。 
随 着 气候 变化 ,我 国 冻 土 大 致 体现 在 :最 大 冻 土 深 
度 下 降 .冻结 日 期 推迟 消融 日 期 提前 \ 冻 土 下 限 上 
升 “"。20 世 纪 70 年 代 至 21 世 纪 初 ,我 国 冻 土 面积 
减少 约 18.6% ,面积 由 2.15x10' km 减少 到 1.75x105 
km 。 按 照 对 未 来 的 天 气 情况 预报 ,未 来 冰冻 圈 各 
组 成 部 分 面积 将 持续 减少 ,而 多 年 冻 土 的 活动 层 厚 
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度 将 会 进一步 增加 ,季节 冻 土 的 冻结 深度 也 将 进 一 
PTE, 


冻 融 指数 是 研究 冻 土 特性 的 重要 参数 ,也 是 判 
定 其 存在 和 长 期 变化 的 一 个 重要 指标 ""。 冻 土 研 
究 中 冻 融 指数 有 助 于 绘制 和 预测 永久 冻 土 分 布 , 佑 
算 土壤 冻结 深度 ,并 计算 活跃 层 厚度 ,冻结 指数 的 
下 降 和 融化 指数 的 上 升 在 一 定 程度 上 说 明 冻 土 环 
境 的 退化 "”。 利 用 季节 冻 土 区 的 冻结 指数 ,可 以 预 
测 季节 冻 融 深度 3; 利用 多 年 冻 土 区 的 冻 融 指数 ， 
不 仅 能 够 绘制 出 冻 土 的 分 布 图 ,还 能 对 多 年 冻 土 的 
活动 层 厚度 进行 预测 “2 ;Peng 等 呈 的 研究 结果 显 
AR ,北半球 高 纬度 和 高 海拔 地 区 的 冻结 指数 下 降幅 
度 最 大 ,冻结 指数 的 下 降 速率 要 比 融 化 指数 的 上 升 
速率 快 得 多 ;Shi 等 "研究 表明 ,从 1901 一 2015 年 ， 
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环北 极地 区 的 冻结 和 融化 指数 分 别 呈 现 降低 、 上 升 
趋势 ,在 1988 年 之 后 ,其 变化 趋势 更 为 明显 ;Liu 等 ” 
基于 地 面 温度 观测 资料 和 再 分 析 产 品 计算 了 雅 鲁 
藏 布 江 上 游 冻 融 指数 ,并 基于 冻 融 指数 模拟 了 冻 土 
面积 的 变化 ,效果 良好 ;Wu 等 5 对 青藏 高 原 冻 融 指 
数 的 变化 规律 进行 了 分 析 ,发 现 多 年 冻 土 区 的 冻结 
指数 比 季节 冻 土 区 的 冻结 指数 降低 速度 要 快 ,而 融 
化 指数 上 升 速 率 更 慢 ;Wang 等 "研究 认为 ,黄河 源 
区 海拔 高 度 对 冻 融 指数 的 空间 分 布 有 一 定 的 影响 ， 
地 表 冻 融 指数 的 变化 比 空气 冻 融 指数 的 变化 要 大 。 

内 蒙古 地 区 的 气象 观测 站 数量 较 少 ,有 些 站 点 
建设 时 间 较 晚 ,总 体 分 布 稀 蚊 ,个 别 站 点 存在 数据 
缺失 问题 ,气候 变化 条 件 下 的 冻 融 指数 变化 研究 较 
为 欠缺 。 仅 有 的 研究 主要 聚焦 在 典型 草原 或 多 年 
冻 土 区 , 且 研 究 对 象 为 空气 冻 融 指数 ”2” ,不 能 完全 
代表 内 蒙古 全 区 的 地 表 冻 融 特 征 ,并 且 研 究 区 地 表 
冻 融 指数 变化 的 驱动 因素 有 竺 探究。 鉴于 此 ,本文 
基于 内 蒙古 气象 站 1980 一 2019 年 日 均 地 表 温 度数 
据 结合 中 国 第 一 代 全 球 陆 面 再 分 析 产 品 (CRA ) 数 
据 以 及 NDVI 数 据 ,对 内 蒙古 近 40 a 地 表 冻 融 指 数 
时 空 变化 特征 及 驱动 因素 进行 分 析 ,该 研究 能 够 为 
下 一 步 估算 内 蒙古 季节 冻 土 冻结 深度 、 多 年 冻 土 活 
动 层 厚度 .模拟 冻 土 面积 变化 等 提供 重要 参数 数 
据 ,以 及 为 内 蒙古 农 牧 业 生产 布局 .草场 生态 评价 
及 冻害 防治 等 方面 提供 一 定 的 科学 依据 。 


100°E 110°E 


120°E 


(a) 气象 站 点 / 冻 土 区 分 布 


a 

j 海拔 Im 。 气象 站 点 
3526 多 年 冻 土 区 

Mm, © stKLK 


0 250km 
| es | 


100°E 110°E 120°E 


注 : 底 图 采用 


50°N 


40°N 


张 吴 琛 等 :内 蒙古 地 表 冻 融 指数 动态 变化 与 驱动 因素 分 析 


ChinaXiv 合 作 期 刊 


1997 


1 数据 与 方法 


1.1 研究 区 概况 

内 蒙古 位 于 126*04' ~ 97° 12'E,37°24' ~53°23'N 
之 间 ( 图 1), 地 处 北部 边疆 ,总 面积 118.3x10'km?, 占 
全 国 总 面积 的 12.3%。 地 貌 以 高 原 为 主 , 兼 有 平 
原 . 沙 地 .丘陵 与 山地 ,大 部 分 地 区 海拔 在 1000 m 
以 上 2 。 气 候 属 温 带 大 陆 性 季风 气候 ,处 于 干旱 与 
半 干 旱 气 候 带 ,年 平均 气温 0~8 %C, 且 由 东北 向 西南 
递增 ;全 区 降水 量 少 日 分布 不 均 , 年 总 降水 量 50~ 
450 mm, 由 东北 向 西南 递减 池 。 气 候 分 异 明显 使 得 
人 研究 区 域 从 东北 向 西南 呈现 出 寒 温 带 针 叶 林 、 落 
叶 闪 时 林 . 草 多 草原 .典型 草原 re Ye Jk AL BK 
分 布下 。 人 研究 区 位 于 环北 极地 区 多 年 冻 土 区 欧 亚 
大 陆 多 年 冻 土 南 界 。 冬 季 严 寒 , 夏 季 高 温 , 且 局 地 
日 夜 温差 大 , 春 冬 季 土 壤 冻 融 现 象 突出 ,夏秋 季 圭 
EPELE”. ZEREZ THEN 
尔 地 区 ,其 余地 区 广泛 分 布 季节 冻 土 , 浅 层 表层 十 
FEA ZAR UK ,夏季 完全 解冻 。 
1.2 数据 来 源 与 处 理 

使 用 的 45 个 气象 台 站 1980—2019 年 日 平均 0 
cm 地 表 温 度 观测 数据 来 自 于 国家 气象 信息 中 心 的 
中 国 地 面 气候 资料 日 值 数 据 集 Chttp://data.cma. 
cn/)。 为 更 好 地 体现 研究 区 的 冻 融 特征 ,尽量 采用 
过 去 40 a 来 不 同 观测 站 的 共同 时 间 数 据 来 求 取 地 表 
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自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2019)1822 号 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 内 蒙古 气象 站 点 / 冻 了 


区 分 布 及 土地 歼 盖 类 型 


Fig. 1 Distribution of meteorological stations/permafrost regions and land cover types in Inner Mongolia 
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冻 融 指数 。 对 缺 测 数据 采用 如 下 方式 内 插 下 : 缺 测 
1 d, 对 其 前 后 各 1 4 数据 的 平均 值 内 搬 ; 缺 测 2 d, ZE 
缺 测 当天 选择 该 日 前 2 d 的 数值 取 平 均 , 缺 测 第 2 d 
选取 该 日 后 2 d 的 数值 取 平均 ;如 果 有 连续 2 4 以 上 
的 缺 测 ,可 以 利用 对 相 邻 站 同一 月 份 的 温度 平均 值 
做 线性 相关 性 检查 ,再 选取 温度 相关 性 最 佳 的 站 点 
用 线性 回归 的 温度 估计 值 内 插 。 采 用 反 距离 插值 法 
(inverse distance weighted) 对 内 蒙古 45 个 气象 站 点 的 
地 表 冻 融 指数 进行 空间 插值 ,插值 分 辨 率 为 34 kmo 

冻 土 类 型 分 区 依据 青藏 高 原 科 学 数据 中 心 
(http://data.tpdc.ac.cn) 提 供 的 1:1000 冻 土 分 区 和 类 型 
地 图 资料 进行 划分 ,通过 对 冻 土 类 型 整合 ,整体 划 
分 为 多 年 冻 土 区 和 季节 冻 土 区 ,用 来 对 地 表 冻 融 指 
数 和 影响 因素 分 区 讨论 。0.1 m 和 0.4 m KERERE 
水 量 、 降 水 量 、 雪 深 数 据 均 源 自我 国 国家 气象 信息 
中 心 的 中 国 第 一 代 全 球 大 气 和 陆 面 再 分 析 产 品 
CRA/Land- 逐 月 陆 面 产品 (http:/data.cma.cn/) ,空间 
分 辨 率 为 34 km, 用 来 分 析 内 蒙古 地 表 冻 融 指数 影 
响 因素 时 空 变化 。 

GIMMS NDVI3g 数 据 由 美国 马里 兰 大 学 GLCF 
(Global Land Cover Facility) 人 研究 组 生产 的 15 d 合 成 
产 品 数据 (https://ecocast.arc.nasa.gov/data/pub/ 
gimms/), MODIS NDVI 数据 为 美国 国家 航天 航空 局 
(NASA) 提 供 的 MODIS 16 d 合 成 产品 数据 (https://mo- 
dis.gsfc.nasa.gov/data/)。 在 GIMMS NDVI 和 MODIS 
NDVI 的 基础 上 ,采用 线性 回归 分 析 , 对 2 种 数据 进 
行 时 间 序 列 的 重建 ,从 而 使 数据 重 采 样 到 同一 个 时 
空 尺 度 上 。 使 用 重合 时 间 的 GIMMS NDVI 和 MO- 
DIS NDVI 进 行 建 模 , 根 据 建 模 结 果 对 GIMMS 数 据 
作 延 长 ,获取 1982 一 2019 年 NDVI 数 据 。 对 于 NDVI 
数据 的 预 处 理 , 主 要 分 为 坐标 变换 .数值 裁剪 ,数据 
转换 .异常 值 删除 .月 最 大 值 获取 等 过 程 , 使 参与 计 
算 的 所 有 NDVI 值 在 0~1 范 围 内 ,根据 研究 区 土地 
覆盖 分 类 图 (图 1) ,结合 研究 区 实际 情况 ,认为 ND- 
VI<0.15 的 地 区 无 植被 覆盖 。 考 虑 到 数据 空间 分 辨 
率 的 一 致 性 ,对 NDVI 数 据 重 采样 为 34 kmo 
1.3 研究 方法 

(1) 地 表 冻 融 指 数 计算 。 基 于 Klene 等 和 
Frauenfeld 等 对 冻结 其 和 融化 期 的 限制 ,界定 冻结 
期 介 于 每 年 7 月 1 日 至 次 年 6 月 30 日 ,融化 期 介 于 
每 年 1 月 1 日 至 12 月 31 日 ,以 最 大 可 能 统计 一 年 中 
寒冷 期 (温度 <0 °C ) 和 温暖 期 (温度 >0 C ) 的 温度 累 


加 值 ,并 确保 负 温 和 正 温 计 算 分 别 在 整个 寒冷 期 和 
温暖 期 内 。 年 冻结 指数 计算 如 下 : 


N, 
FI= >it T,<0°C 
i=l 


(1) 


N; 
= 7T>0 C (2) 
= 


SUE FLAN TI PG RG HS ETEC d); T, 
和 了 分 别 为 负 值 逐 日 地 表 温 度 和 正 值 逐 日 地 表 温 
BE, ALAC; N, AlN, 为 年 内 地 表 温 度 小 于 和 大 
于 0% 的 天 数 ,单位 为 d。 

(2) 趋势 分 析 。 利 用 一 元 线性 回归 方法 计算 冻 
融 指数 和 影响 因素 的 年 际 变化 趋势 ,并 通过 最 小 二 
乘法 拟 合 的 斜率 反映 , 当 和 斜率 大 于 0 时 ,说明 有 增 大 
的 倾向 , 当 斜 率 小 于 0 时 ,说 明 有 下 降 的 倾向 。 其 变 
化 趋势 的 计算 公式 为 : 

nix X-ið, 


i=l i=l i=l 


Slope = 


式 中 :Slope 表示 地 表 冻 融 指 数 或 影响 因素 的 变化 趋 
势 ;n 为 时 间 序 列 长 度 ; X, HWA RAX T i FIF 
均值 。 

(3) 相关 分 析 法 。 通 过 计算 地 表 冻 融 指数 与 对 
应 的 影响 因素 相关 系数 ,表示 其 与 地 表 冻 融 指数 的 
相关 性 。 计 算 公 式 如 下 : 


G0) 
r= a (4) 


xy 
y n 


Sena So- 
式 中 :n 表 示人 研究 时 间 段 的 累计 年 数 ; r, 表示 x 和 y 
的 相关 系数 ; x, 和 yy, 分 别 表示 第 i 年 x 和 y 的 值 ; x 
和 了 表示 2 个 要 素 样本 值 的 平均 值 ; 7, 在 -1~1 之 
间 , 随 着 绝对 值 的 增 大 ,地表 冻 融 指数 与 影响 因素 
的 变化 关系 也 相应 较 好 ;反之 越 差 , 当 7，>0 表 示 
相关 性 为 正 , r,。<0, 表 示 相 关 性 为 负 ”。 最 后 采用 
了 检验 法 进行 显著 性 检验 ,将 地 表 冻 融 指 数 与 影响 
因素 的 相关 程度 划分 为 显著 相关 (0.01< P<0.05) 
和 不 显著 相关 (P 宇 0.05)。 

(4) 灰色 关联 分 析 。 采 用 灰色 关联 分 析 法 , 研 
究 影响 地 表 冻 融 指数 的 主要 因子 及 重要 性 。 灰 色 
关联 度 可 以 测量 2 个 系统 之 间 的 关系 ,并 确定 相关 
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子 中 的 主 控 因 子 ”。 地 表 冻 融 指数 的 时 间 序 列 

作为 目标 变量 ,将 0.1 m 土壤 含水 量 、0.4 m 土壤 含 水 

量 \ 降 水 量 等 影响 因素 作为 影响 变量 。 灰 色 关 联系 

数 表达 为 : 

内 =- min, min,| ylk) — x(k) | +p max, max,| y(k)- «,(h)| 
|y) - x;(k)| +p max, max,| y(k) — x,(h)| 


(5) 
G.= 1 (6) 
k=1 


式 中 : é, (有 ) 为 影响 变量 与 目标 变量 的 灰色 关联 
系数 ;po 为 分 辩 系 数 ,参考 以 往 研 究 ,在 本 文中 p 取 
0.5 2 ; min, miny) -x| 为 二 级 最 小 差 ; 
p max, max,| ylk) — x,(k)| H ZRARKZŽ ; | yk) -x| A 
各 比较 数列 x, 曲线 上 的 每 个 点 与 参考 数列 x, 曲线 
上 每 个 点 差 的 绝对 值 。 6G 为 比较 数列 x, 对 参考 数 
列 xu 的 灰色 关联 度 ,用 以 表达 二 者 之 间 的 联系 紧密 
程度 ,该 值 越 接近 1 ,说 明 相关 性 越 好 。 灰 色 关 联 度 
最 大 的 影响 变量 为 目标 变量 的 主导 影响 因素 , 即 主 
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2 结果 与 分 析 


21 内 蒙古 地 表 冻 融 指 数 时 空 变 化 

2.1.1 地 表 冻 融 指数 的 时 间 变 化 特征 ”通过 对 内 蒙 
古 地 区 1980 一 2019 年 的 45 个 观测 站 SFI( 地 表 冻 结 
指数 ) 和 STI( 地 表 融 化 指数 ) 进 行 分 析 ,得 到 人 研究 区 
40 a 整体 和 不 同 冻 土 区 SET 与 STI 的 时 间 变 化 特 
征 。 如 图 2a 和 图 2b 所 示 , 整 体 上 SFI 近 40 a I i. 
著 下 降 趋 热 ,多 年 冻 土 区 SFI 下 降 速 率 高 于 季节 冻 
土 区 ,内 蒙古 SFI 多 年 变化 范围 介 于 956.1~ 
1848.3 "Cd, 最 低 值 和 最 高 值 分 别 发 生 于 2016 年 和 
1984 年 ,年 均值 为 1457.0 ‘Cd, 变化 率 为 -156.4 C- 
d.(10a)-(P<0.01) ,表明 土 壤 的 季节 冻结 潜力 有 连 
续 降 低 的 趋势 ”。 进 一 步 分 析 各 站 点 SET 变化 特 
征 , 近 40a 所 有 站 点 的 SFEI 均 呈 下 降 趋 势 , 但 各 站 
SFI 变 化 特征 差异 也 较 大 ,变化 最 小 值 和 最 大 值 分 别 
出 现 于 阿拉 善 盟 阿 拉 善 右 旗 站 点 [-20.7%.d (10a)”] 
和 兴安 盟 阿尔 山 站 点 [-711.1 Cd*(10a)-'] 。 如 图 
2c 和 图 2d 所 示 , 整 体 上 STI 近 40 a 呈现 显著 上 升 趋 
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注 : 虚 线 代 表 相 应 的 线性 趋势 。 下 同 。 
图 2 1980 一 2019 年 内 蒙古 整体 和 不 同 冻 土 区 地 表 冻 融 指数 年 际 变化 趋势 
Fig.2 Interannual variation trend of surface freezing-thawing indices in Inner Mongolia and different permafrost 


regions from 1980 to 2019 


202301.00127v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


2000 干 F 究 39 卷 


El 
3 
È 


势 ,多 年 冻 土 区 STI 上 升 速率 高 于 季节 冻 土 区 ,内 蒙 ”区 位 于 阿拉 善 西 北部 地 区 , STREKT IEE D 
古 STI 多 年 变化 范围 介 于 3717.6~4442.3 Cd, 最 低 。 尔 中 部 地 区 。 内 蒙古 STI 年 际 变 化 率 空间 分 布 特征 
值 和 最 高 值 分 别 发 生 于 1984 年 和 2007 年 ,年 均值 见 图 3d,STI 上 升 速率 较 高 的 地 区 (Slope > 150 Cd: 
为 4040.9 Cd, 上 升 速 率 为 152.4 Cd*(10a)"(P< (10a)) 主 要 分 布 在 研究 区 中 部 .中 东部 和 东北 部 
0.01), 说 明 土 壤 融 化 期 间 最 大 融化 深度 增加 中。 ”分 地 区 ,可 能 是 因为 STT 上 升 速率 较 高 的 地 区 受 人 
1980 一 2019 年 间 , 所 有 站 点 的 STI 均 呈 上 升 趋势 ,各 ”类 活动 等 因素 的 影响 ,其 次 ,呼伦贝尔 高 原 大 部 分 
站 STI 变化 特征 存在 差异 ,变化 特征 最 小 值 和 最 大 地 区 的 冻 融 指数 变化 率 最 大 ,这 也 证 实 了 该 区 域 
值 分 别 出 现 于 赤峰 市 赤峰 站 点 [43.1 "Cd*(10a)"] ”( 原 岛 状 冻 土 区 ) 目 前 处 于 季节 性 冻 土 的 状态 。 
和 赤峰 市 巴林 左旋 站 点 [221.1%.d.(10a)]。 2.2 内 蒙古 地 表 冻 融 指 数 影 响 因素 时 空 变 化 
2.1.2 地 表 冻 融 指数 空间 分 布 与 变化 ”内 蒙古 SFI 22.1 内 蒙古 地 表 冻 融 指 数 影响 因素 时 间 变 化 “如 
多 年 平均 值 的 空间 分 布 特征 如 图 3a 所 示 , 受 经 度 、 图 4 所 示 , 近 40 a 研 究 区 0.1 m 和 0.4m 土壤 含水 量 
纬度 和 海拔 的 综合 影响 ,SFI 均 值 总 体 呈 现 西 低 东 ， 年 均值 分 别 介 于 0.1401~0.1673 m° -m° 和 0.1388~ 
高 .由 西南 向 东北 逐渐 增加 的 特征 ,其 中 SFI 受 纬度 0.1779 mm 之 间 ,0.4m 土 壤 含水 量 略 大 于 0.1 m 
的 影响 最 大 ( 表 1) ,相关 系数 达到 0.92(P< 0.01) ,每 土壤 含水 量 , 不 同 深度 土壤 含水 量 的 总 体 变化 趋势 
EFF 1°, SFI SG HN 165.2 %.d, 该 分 布 特征 同 环北 ” 均 显 示 出 向 上 态势 ,平均 年 际 变化 率 为 0.001 mem? 
极 冻 融 指 数 分 布 变 化 特征 一 致 "*" 。SFI 高 值 区 (10a) ,而 0.1m 和 0.4m 土 壤 含水 量 并 不 是 呈现 一 
位 于 呼伦贝尔 北部 地 区 ,SFI 低 值 区 位 于 阿拉 善 。 致 性 的 向 上 态势 ,而 是 随 气候 波动 经 历 了 “上 升 - 下 
西部 地 区 ,与 冻 土 区 分 布 状况 类 似 ,多 年 冻 土 区 降 - 上 升 "3 个 波动 性 变化 阶段 ,是 0.1m 和 0.4m 深 
SFI 普 遍 大 于 季节 冻 土 区 。 内 蒙古 SFI 年 际 变 化 。 上 度 的 土壤 含水 量 年 际 变 化 趋势 大 体 相同 :1980 一 
率 空间 分 布 特征 见 图 3b,SFI 年 际 变化 率 的 变化 1998 年 为 波动 上 升 阶段 ,1998 一 2007 年 呈 波 动 下 降 
范围 介 于 -10.3~-710.6%.d. (10a), SIRER ”趋势 ,2007 一 2019 年 又 呈 波 动 上 升 趋势 ;降水 量 以 
在 空间 上 自 西南 向 东北 呈 增 加 趋势 ,表明 内 蒙古 10.0 mme (10a) 的 速率 呈 不 显著 增加 趋势 (P> 
东北 地 区 SFI 下 降 速率 较 西 南 地 区 更 快 。 0.05), 年 均值 介 于 240.1~415.9 mm 之 间 ,降水 量 与 
内 蒙古 STI 多 年 平均 值 的 空间 分 布 特征 见 图 “0.1m 和 0.4m 土 壤 含 水 量 的 年 际 波动 变化 趋势 极为 
3c,STI 均 值 的 空间 分 布 特征 呈现 西部 高 ,东部 低 ， “相似 ,同样 经 历 了 “ 上升- 下降- 上 升 ?3 个 波动 变化 
由 西南 向 东北 逐渐 减 小 的 趋势 ,多 年 冻 土 区 的 STI ”阶段 ,1980 一 1998 年 .2007 一 2019 年 属于 降水 偏 丰 
普遍 小 于 季节 冻 土 区 ,在 呼和浩特 、 赤 峰 等 地 出 现 ”的 年 代 ;1998 一 2007 年 属于 降水 偏 村 的 年 代 ; 研 究 
了 点 状 高 值 区 ,这 与 近 几 年 农业 开发 .城市 扩张 . 工 区 雪 深 以 11.4~32.8 mm 的 起 伏 趋势 为 主 ,多 年 积 雪 
业 发 展 等 社会 经 济 活动 密切 相关 。 同 时 ,上 述 地 区 ”厚度 为 22.1 mm ,年 平均 雪 深 以 -0.82 mm: (10a) HS 
处 于 生态 环境 较为 脆弱 ,人 为 活动 对 环境 的 影响 更 。 速率 呈 微 弱 减 少 趋势 (P > 0.05); 基于 1982—2019 年 
加 敏感 。STI 同 纬度 的 相关 系数 最 高 ( 表 1), 为 -0.73 ”内 蒙古 NDVI 数 据 , 发 现 内 蒙古 NDVI 总 体 以 0.007: 
(P<0.01), 每 上 升 1° STI 减 少 141.1 Cd,STI 高 值 (10a)"' 的 速率 明显 著 增加 趋 热 (P <0.05), 年 均值 介 
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图 3 1980 一 2019 年 内 蒙古 地 表 冻 融 指 数 均值 和 变化 趋势 


Fig. 3 Mean value and variation trend of surface freezing-thawing index in Inner Mongolia from 1980 to 2019 
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R1 地 表 冻 融 指数 与 地 理 因 子 的 相关 系数 
Tab. 1 Correlation coefficient between surface freezing- 


thawing index and geographical factors 


因子 年 均 SFI 年 均 STI 
经 度 0.63” -0.62" 
纬度 0.92" -0.73" 
观测 场 拔 海 高 度 /m -0.38" 0.12 


N 


注 :* 表 示 通 过 P<0.05 显著 水 平 检 验 ;#*# 表 示 通 过 P<0.01 显著 水 平 


于 0.41~0.47 之 间 。 

2.2.2 内 蒙古 地 表 冻 融 指 数 影响 因素 空间 分 布 与 变 
化 ”如 图 5 所 示 , 研 究 区 内 各 影响 因素 空间 分 布 整 
体 上 存在 明显 的 纬度 地 带 性 和 区 域 性 规律 , 且 均 呈 
现 出 由 西南 向 东北 递减 的 发 展 趋势 ,与 SET 空间 分 
布 特征 一 致 ,与 STI 空 间 分 布 特征 相反 。 但 上 述 要 
素 在 部 分 地 区 空间 分 布 和 空间 变化 存在 一 定 差异 : 
0.1 m 和 0.4 m 土 壤 含 水 量 呈 现 出 东北 部 和 中 西部 土 
二 水 分 相对 较 高 .西部 相对 较 低 的 空间 格局 ,可 能 
受到 不 同 土地 覆盖 类 型 . 积 雪 分 布 及 人 为 活动 等 多 
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重 因素 的 共同 作用 ;由 于 大 兴安 岭 .阴山 、 贺 兰 山 等 
地 势 因素 的 阻碍 ,东南 地 区 的 水 汽 抬升 ,容易 造成 
降雨 ,西北 山区 的 水 汽 含量 偏 低 ,因此 东南 地 区 的 
降水 量 比 西北 地 区 大 ; 雪 深 与 降水 量 分 布 特征 基本 
一 致 ,但 与 温度 相反 ,研究 区 东北 部 降水 量 较 多 AF 
均 地 表 温 度 较 低 且 森林 有 覆盖 率 高 ,有 利于 积 雪 的 形 
成 和 积累 ,中 部 和 西部 地 区 降水 较 少 ,平均 地 表 温 
度 较 高 ,植被 覆盖 率 较 低 ,不 利于 积 雪 的 形成 与 累 
R ,使 得 雪 次 较 浅 且 分 布 均匀 ;内 蒙古 NDVI 由 西南 
向 东北 呈 增 加 趋势 ,部 分 地 区 NDVI 存 在 区 域 性 特 
征 和 空间 异 质 性 较 明 显 , 可 能 与 灌溉 和 开明 农 田 等 
人 类 活动 有 关 。 

如 图 6 所 示 , XE 40 a 来 内 蒙古 大 部 分 地 区 0.1 m 
和 0.4 m 土 壤 含 水 量 年 际 变化 呈 上 升 态势 ,0.1 m 和 
0.4 m 土 壤 含 水 量 呈 减少 态势 的 区 域 分 别 仅 占 研 究 
区 整体 的 19.9% 和 36.1%, 且 大 部 分 集中 在 研究 区 西 
北部 和 东北 部 的 多 年 冻 土 区 ;内 蒙古 中 部 和 西部 的 
年 降水 量 有 所 增加 ,特别 是 中 偏 南 地 区 降水 量 增加 
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图 4 1980 一 2019 年 影响 因素 年 际 变化 趋势 


Fig.4 Interannual variation trend of impact factors from 1980 to 2019 


202301.00127v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


2002 干旱 区 研究 39 4 
100*E 110°E 120°E 100°E 110°E 120°E 100°E 110°E 120°E 
(a) 0.1 m 土 壤 含水 量 A (b) 0.4 m 土 壤 含水 量 z 
Zz 7 o y o Z. ° 
E 上 土壤 含水 量 (mm B E 土壤 合 水 量 /mm 3 È i 
= 0.286 = 0.301 
0.044 0.071 
Zz z z 
z S Z 2 2 & 
S FS +S + 
T <+ + 
100°E 110°E 120°E 100°E 110°E 120°E 100°E 110°E 120°E 
100°E 110°E 120°E 100°E 110°E 120°E 
(d) 雪 深 ? (e) NDVI 这 
Z j: >o Z > 
b 积 雪 深度 /mm a è NDVI 值 aS 
137.1 0.92 
= = 号 无 植被 区 
0.3 0.15 
rA rA 
£ $Z 2 
时 S 
100°E 110°E 120°E 100°E 110°E 120°E 


图 5 1980 一 2019 年 地 表 冻 融 指 数 影 响 因素 多 年 平均 


Fig. 5 Annual average of surface freezing-thawing index influencing factors from 1980 to 2019 
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图 6 1980 一 2019 年 地 表 冻 融 指数 影响 因素 年 际 变化 


Fig.6 Interannual variation of surface freezing-thawing index influencing factors from 1980 to 2019 


较 多 ,东部 降水 量 则 有 所 减少 ,特别 是 在 呼伦贝尔 
西北 部 .赤峰 中 部 以 及 通辽 中 部 地 区 降水 量 减 少 情 
况 较为 突出 ;研究 区 38.5% 的 区 域 雪 深 呈 下降 态势 ， 


FL 10.8% 的 区 域 雪 深 平 均 下 降 速 度 在 0.5 mm: 
(10a)"' 以 上 ,主要 分 布 在 呼伦贝尔 东北 部 ,并 位 于 
多 年 冻 土 区 ;61.5% 的 区 域 雪 深 呈 现 增 加 趋势 ,主要 
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分 布 于 西南 部 的 沙漠 及 中 部 和 东部 的 部 分 草原 ,多 
位 于 季节 冻 土 区 , 且 以 上 区 域 的 雪 深 平均 增加 速率 
仅 为 0.4 mm (10a)" ,这 一 变化 特征 可 能 与 研究 区 
近年 来 温度 和 降水 变化 有 关 ; 研究 时 段 内 蒙古 有 
75.2% 区 域 的 NDVI 呈 现 增 长 的 趋势 , 仅 24.8 色 的 区 
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表 3 SFI 与 影响 因素 相关 系数 面积 占 比 
Tab.3 Area ratio of correlation coefficient between SFI 


and influencing factors 


相关 性 ”0.1m 土 壤 0.4 m5 


含水 量 含水 量 


降水 量 S% NDVI 


域 呈现 下 降 的 趋势 ,下 降 区 域 以 锡林郭勒 盟 大 部 、 
赤峰 市 北部 .通辽 市 北部 .兴安 盟 南 部 .呼伦贝尔 市 
林 区 为 主 ,这 一 变化 特征 可 能 与 上 述 地 区 近年 来 温 
度 升 高 .降水 减少 及 人 类 活动 有 关 。 
2.3 地 表 冻 融 指 数 与 影响 因素 相关 分 析 

地 表 冻 融 指数 的 变化 受 多 种 因素 的 影响 ,具有 
较 强 的 随机 性 , 适 于 进行 灰色 关联 度 分 析 。 如 表 2 
所 示 ,内 蒙古 地 表 冻 融 指 数 年 均值 与 0.1 m 土壤 会 
水 量 、0.4 m 土 壤 舍 水 量 、 降 水 量 、 雪 深 、NDVI 的 灰色 


显著 正 相关 1.1% 8.9% 0.0% 5.3% 3.3% 
不 显著 正 相 关 16.0% 24.0% 11.4% 43.2% 20.8% 
不 显著 负 相 关 35.0% 29.9% 55.9% 46.3% 48.3% 
显著 负 相关 47.9% 37.2% 32.7% 5.2% 24.6% 


影响 因素 相关 性 差异 的 原因 可 能 与 影响 因素 的 时 
空 变 化 差异 有 关 。 其 中 ,0.1m 土壤 含水 量 与 SFI 显 
著 负 相关 (P<0.05) 的 区 域 大 部 分 在 季节 冻 土 区 的 中 
部 和 东南 部 ,显著 正 相 关 (P<0.05) 区 域 集中 分 布 在 


关联 度 系数 缘 大 于 0.65 ,说 明 地 表 冻 融 指数 变化 与 
影响 因素 存在 较 大 的 相关 性 。SFI 与 影响 因素 的 
灰色 关联 序 为 :0.4 m 土 壤 含 水 量 > 0.1 mm 土壤 含水 
量 >NDVI> 降水 量 > 雪 深 , 说 明了 0.4m 土 壤 含 水 
量 对 SFI 变 化 的 影响 较 大 ,是 影响 内 蒙古 整体 SFI 变 
化 的 主 控 因素 ,0.1 m 土壤 含水 量 .NDVI、 降 水 量 次 
之 ,最 后 是 雪 深 ;STI 与 影响 因素 的 灰色 关联 序 为 : 
NDVI> 0.1 m 土壤 含水 量 > 0.4m 土壤 含水 量 > 降 
水 量 > 雪 深 ,表明 NDVI 对 STI 变 化 的 影响 较 大 ,是 
影响 内 蒙古 整体 STI 变 化 的 主 控 因素 ,0.1 m 土壤 含 
水 量 、0.4 m 土 壤 售 水量、 降水 量 次 之 ,最 后 是 雪 深 。 


表 2 地 表 冻 融 指数 与 影响 因素 灰色 关联 系 分 析 
Tab.2 Grey correlation analysis between surface 
freezing-thawing index and influencing factors 

地 表 冻 融 指 数 关联 度 结果 
SFI 关 联 度 ”SFI 关 联 序 STI 关 联 度 STI 关 联 序 


影响 因素 


0.1 m 土 壤 湿 度 0.76 2 0.859 2 
0.4 m 土 壤 湿 度 0.761 1 0.833 3 
降水 量 0.685 4 0.739 4 
雪 深 0.681 5 0.671 5 
NDVI 0.754 3 0.864 1 


为 探讨 地 表 冻 融 指数 对 影响 因素 变化 的 动态 
响应 ,基于 像 元 尺度 逐 像 元 计算 了 内 蒙古 地 表 冻 融 
间 数 与 影响 因素 的 相关 关系 。SFI 与 各 影响 因素 
(除雪 深 外 ) 相 关 关 系 有 超过 65% 的 地 区 呈 负 相关 
( 表 3) ,如 图 7 所 示 , 从 空间 分 布 来 看 ,在 多 年 冻 土 区 
大 部 分 地 区 呈现 正 相 关 关 系 , 但 在 季节 冻 土 区 大 部 
分 地 区 却 呈 现 负 相关 关系 ,造成 不 同 冻 土 区 SFI 与 


多 年 冻 土 区 北部 小 部 分 地 区 , 0.4 mm 土壤 含水 量 与 
SFI 显 著 负 相关 (P<0.05) 的 地 区 在 季节 冻 土 区 的 南 
部 和 东南 部 均 有 分 布 , 显 著 正 相关 (P<0.05) 的 地 区 
则 在 多 年 冻 土 区 北部 呈现 片 状 分 布 ;降水 量 与 SFI 
显著 负 相 关 (P<0.05) 的 地 区 主要 在 鄂尔多斯 市 以 
东 、 呼 和 浩特 市 和 锡林郭勒 盟 以 南 ,呈正 相关 关系 
的 区 域 主要 在 多 年 冻 土 区 北部 和 阿拉 善 中 西部 ; 
NDVI 与 SFI 显 著 负 相关 (P<0.05) 的 区 域 主要 在 鄂 尔 
多 斯 东部 .呼和浩特 市 南部 和 多 年 冻 土 区 部 分 林 
地 ,显著 正 相 关 (P<0.05) 区 域 主要 在 多 年 冻 土 区 部 
分 林地 和 湿地 地 区 ; 雪 深 与 SFI 显 著 负 相关 (P<0.05) 
的 区 域 主要 在 阿拉 善 西南 部 和 季节 冻 土 区 中 部 地 
区 零星 分 布 , 显 著 正 相关 (P<0.05) 区 域 集中 在 多 年 
冻 土 区 北部 和 阿拉 善 西北 部 小 片 地 区 。 

STI 与 影响 因素 的 相关 关系 以 负 相 关 为 主 ( 表 
4)。 如 图 8 所 示 , 虽 然 STI 与 各 影响 因素 在 多 年 冻 土 
区 及 大 部 分 地 区 呈现 负 相 关 关 系 , 但 在 季节 冻 土 区 
部 分 地 区 却 呈 现 正 相关 关系 ,造成 不 同 冻 土 区 STI 
与 影响 因素 相关 性 差异 的 原因 可 能 与 影响 因素 的 
时 空 变化 差异 有 关 。 其 中 ,0.1 m 土壤 含水 量 与 STI 
显著 负 相 关 (P<0.05) 区 域 集中 在 多 年 冻 土 区 东部 和 
多 年 冻 土 区 周围 的 盟 市 ,显著 正 相 关 (P<0.05) 的 区 
域 在 黄河 流域 内 蒙古 地 区 有 所 分 布 ;0.4 m 土壤 含水 
量 与 STI 显 著 负 相关 (P<0.05) 区 域 主要 在 多 年 冻 土 
区 北部 和 阿拉 善 西部 零星 地 区 ,显著 正 相 关 (P< 
0.05) 的 区 域 主要 在 阿拉 善 东部 E ET TH Aaa A 
鄂尔多斯 市 周边 ;降水 量 与 STI 显 著 负 相关 (P<0.05) 
区 域 主要 在 锡林郭勒 盟 东 部 .兴安 盟 东部 及 多 年 冻 
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7 SFI 与 影响 因素 相关 性 7 检验 


Fig. 7 Correlation T test between SFI and influencing factors 


#4 STI 与 影响 因素 相关 系数 面积 占 比 
Tab.4 Area ratio of correlation coefficient between STI 
and influencing factors 


STI 

ere sae ea pee 降水 量 ， 雪 深 NDVI 
显著 正 相 关 8.1% 12.1% 0.0% 2.0% 18.1% 
不 显著 正 相 关 36.9% 24.5% 34.1% 35.8% 33.8% 
不 显著 负 相 关 31.5% 29.0% 53.8% 54.5% 32.2% 
显著 负 相 关 23.5% 34.4% 12.1% 7.7% 12.9% 


土 区 东北 部 ,呈正 相关 关系 的 区 域 主要 在 季节 冻 土 
区 中 部 和 西部 ; 雪 深 与 STI 显 著 负 相关 (P<0.05) 区 域 
主要 在 锡林郭勒 盟 西 北部 、 赤 峰 市 北部 及 多 年 冻 土 
区 东北 部 ,显著 正 相关 (P<0.05) 的 区 域 集 中 在 巴 彦 
淖尔 市 中 部 。NDVI 与 SFI 显 著 负 相关 (P<0.05) 的 区 
域 主要 在 赤峰 市 北部 、 兴 安 盟 、 多 年 冻 土 区 部 分 林 
地 和 锡林郭勒 盟 车 漠 草 原 地 区 ,显著 正 相 关 (P< 
0.05) 的 区 域 主 要 在 黄河 流域 内 蒙古 地 区 .通辽 市 和 
多 年 冻 土 区 北部 。 


3 讨论 


气候 变 暖 背景 下 ,内 蒙古 1980 一 2019 年 间 地 表 


冻结 和 融化 指数 分 别 显著 上 升 和 下 降 ,多 年 冻 土 区 
地 表 冻 融 指 数 变 化 速率 较 季 节 浆 土 区 更 快 。 人 研究 
结果 基本 与 Peng 等 Wu 等 EARP ILE 
北半球 .青藏 高 原 .东北 地 区 的 结论 相 吻 合 ,说 明 在 
内 蒙古 整体 冻 土 环境 退化 明显 ,多 年 冻 土 区 冬季 温 
度 变 暖 趋势 更 为 显著 ,其 冻 土 退化 更 严重 ,该 区 地 
表 冻 融 指数 的 显著 变化 可 能 导致 多 年 冻 土 区 面积 
缩减 ,冻结 深度 减少 ,活动 层 厚度 增加 等 一 系列 冻 
土 退 化 现象 的 发 生 ” ;季节 冻 土 区 的 地 表 冻 融 指数 
变化 速率 同样 不 容 忽视 ,该 区 是 内 蒙古 农 牧 业 生产 
主要 布局 地 区 , 冻 融 指数 的 变化 同 农 作物 的 种 植 ， 
畜牧 业 的 发 展 息息相关 。 和 气候 变 暖 虽然 是 造成 该 
区 多 年 冻 土 区 冻 土 退化 显著 的 根本 性 因素 ,但 是 多 
种 人 为 经 济 活动 对 研究 区 多 年 冻 土 退化 起 到 了 加 
速 促 进 作用 ”。 近 40 a 人 研究 区 表层 土壤 含水 量 、 降 
水 量 的 增加 趋势 表明 ,内 蒙古 地 区 的 冻 土 环境 正 向 
暖 湿 化 发 展 ,与 马 梓 策 等 ”得 出 的 研究 结论 相 一 
致 。 但 是 基于 本 文 研究 发 现 ,不 同 冻 土 区 的 发 展 情 
况 存在 较 大 差异 ,多 年 冻 土 区 正 向 暧 干 化 发 展 ,与 
研究 区 整体 和 季节 冻 土 区 大 部 分 地 区 变化 趋势 相 
反 。 地面 温度 升 高 降水 减少 ,潜在 蒸发 量 增加 是 
导致 土壤 含水 量 下降 的 主要 因素 ”。 地 表 冻 融 指 
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图 8 STI 与 影响 因素 相关 性 7 检验 


Fig. 8 Correlation T test between STI and influencing factors 


数 是 根据 地 表 温 度 计算 所 得 ,与 气温 和 地 表 温 度 有 
密切 的 联系 ,但 表层 土壤 含水 量 KEk E, SRA 
NDVI 同 样 对 冻 融 指数 的 变化 起 着 至 关 重 要 的 作 
用 ,总 体 上 ,表层 土壤 含水 量 和 NDVI 对 地 表 冻 融 指 
数 的 变化 起 主导 作用 。 此 外 , 冻 土 环境 退化 对 生物 
多 样 性 的 作用 是 负面 的 。 温 度 增高 .雨量 减少 的 情 
况 对 草原 植物 生长 发 育 不 利 ,对 森林 ,湿地 植物 有 
促进 作用 3; 季节 冻 土地 带 冬季 增 温 明 显 ,农作物 
冻害 损失 减少 ,从 而 促进 了 粮食 作物 单产 的 增加 和 
质量 提高 ,利于 家 冀 健 康 过 冬 度 春 ,同时 也 会 使 病 虫 
害 及 虫 卵 易于 过 冬 ,容易 造成 病虫害 的 防治 困难 ”。 
20 世 纪 80 年 代 前 研究 区 气象 数据 缺 测 较 多 ,本 
研究 只 采用 了 1980 一 2019 年 的 气象 资料 ,下 一 步 还 
需 依 据 更 长 时 间 序 列 和 更 多 的 站 点 数据 开展 研 
究 。 在 对 地 表 冻 融 指 数 处理 后 ,需要 根据 冻 融 指 
数 ,进一步 确定 研究 区 多 年 冻 土 空间 分 布 范围 和 活 
动 层 厚度 及 季节 冻 土 冻结 深度 的 时 空 变化 特征 ,对 
内 蒙古 地 区 冻 土 对 气候 变化 的 响应 作出 较为 全 面 、 
精确 的 预测 。 目 前 ,由 于 中 国 全 球 陆 面 再 分 析 产 品 
(CRA/Land) 数 据 在 冻 土 研究 中 仍 具 有 空间 分 辨 率 略 
小 的 缺陷 ,在 今后 冻 土 变化 研究 中 , 拟 将 与 空间 分 
辩 率 更 高 的 再 分 析 资 料及 实际 观测 资料 相 结合 ,以 


尽 可 能 降低 研究 结论 的 不 确定 性 。 
4 结论 


本 文 基 于 内 蒙古 气象 站 1980 一 2019 年 日 均 地 
表 温 度数 据 结合 中 国 第 一 代 全 球 陆 面 青 分 析 产 品 
(CRA) SHH LA NDVI AC He (EB) GIS 空间 分 析 功 
能 ,通过 RStudio 编程 处 理 分 析 了 内 蒙古 近 40 a 地 表 
冻 融 指数 时 空 变化 格局 及 关键 影响 要 素 ,主要 结论 
如 下 : 

(1) 1980—2019 年间, 内 蒙古 SFI 与 STI 分 别 呈 
现 明显 减少 与 增加 的 态势 ,纬度 是 制约 内 蒙古 地 表 
冻 融 指数 空间 分 布 的 重要 因素 ,与 SET 呈现 极 明 显 
正 相 关 的 关系 , 则 与 STT 相 反 。 多 年 冻 土 区 地 表 冻 
融 指 数 变化 速率 较 季 节 冻 土 区 更 为 剧烈 ,对 气候 变 
暧 的 响应 更 为 敏感 。 

(2) 近 40 a 人 研究 区 0.1 m 和 0.4 m EREKE, 
降水 量 .NDVI 均 呈 增 加 趋势 , 仪 雪 深 呈 减少 趋势 。 
其 中 不 同 深度 土壤 含水 量 平均 年 际 变化 率 为 0.001 
mm (10a); KEK LI 10.0 mm: (10a)! AYR SE 
不 显著 增加 趋势 ;NDVI 总 体 以 0.007* (10a) HIRE 
旺 显著 增加 趋势 ;年 平均 雪 深 以 -0. 82 mm (10a)! 
的 速率 呈 微 弱 减 少 趋势 。 
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(3) 多 年 冻 土 区 0.1 m 和 0.4 m 土 壤 含 水 量 、 降 
水 量 .NDVI 及 雪 深 的 多 年 均值 普遍 大 于 季 市 冻 土 
区 ,同样 影响 因素 在 年 际 变化 上 也 表现 出 不 同 的 空 
间 差 异性 ,多 年 冻 土 区 呈 暧 干 化 发 展 趋 势 ,季节 冻 
LK ERA 

(4) 地 表 冻 融 指 数 与 影响 因素 以 负 相 关 关 系 为 
主 ,SFI 与 影响 因素 在 多 年 冻 土 区 大 部 呈正 相关 关 
系 ,在 季 贡 冻 土 区 大 部 呈 负 相关 关系 ,STI 反 之 。 内 
蒙古 地 表 冻 融 指 数 变化 受 影响 因素 共同 驱动 , 而 
0.4m 土 壤 含 水 量 是 影响 SFI 变 化 的 主导 因素 ,NDVI 
是 影响 STI 变 化 的 主导 因素 。 
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Dynamic changes and driving factors of the surface freeze-thaw 
index in Inner Mongolia 
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Abstract: Using the trend analysis, correlation analysis and gray correlation, the spatial and temporal variation 
characteristics and driving factors of the surface freezing index (SFI) and surface thawing index (STI) in Inner 
Mongolia over the past 40 years were analyzed based on the daily average surface temperature data of 45 
meteorological stations in Inner Mongolia from 1980 to 2019, combined with China’ s first- generation global 
land surface reanalysis product (CRA) data and the NDVI data. The study shows that the following: (1) the 
spatial distribution of annual mean values of the SFI had an increasing pattern from southwest to northeast, while 
the STI had the opposite pattern, and latitude was the key factor affecting the spatial distribution of the SFI and 
STI. The SFI and STI showed significant decreasing and increasing trends during the study period, with multi- 
year variations ranging from 956.1 to 1848.3 C -d and 3717.6 to 4442.3 °C-d, respectively, and rates of 
change of —156.4 C-d-(10a)"' and 152.4 C-d-(10a)"'. Compared with the seasonal permafrost zone, the 
freeze-thaw index in the multi-year permafrost zone was more sensitive to climate warming. (2) The soil 
SW content, precipitation, and NDVI in the study area showed an increasing trend, and the snow depth 
showed a decreasing trend over the last 40 years. However, the interannual variation showed different spatial 
variability, with a warm- dry trend in the multi- year permafrost area and a warm- wet trend in the seasonal 
permafrost area. The SFI, STI, and influencing factors are mainly negatively correlated. The SFI and the 
influencing factors were mostly positively correlated in the multi- year permafrost area and mostly negatively 
correlated in the seasonal permafrost area. The STI was the opposite of the SFI. The change of the SFI and STI in 
Inner Mongolia was driven by a combination of influencing factors, while 0.4 m of soil water content was the 
dominant factor affecting the change of the SFI, and NDVI was the dominant factor affecting the change of the 
STI. The results provide valuable references for permafrost degradation and the production of agriculture and 
livestock in Inner Mongolia. 

Keywords: Inner Mongolia; surface freeze- thaw index; climate change; spatio- temporal characteristics; 


driving factors 


